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FOREWORD 


This  work  was  performed  under  US  Army  Natick  Laboratories 
Contract  No,  DAAG  17-72-C-0030  during  the  period  15  November 
1971  and  30  June  1973,  The  project  Number  was  1F162203AA33 
and  the  Task  Number  was  04  entitled  "Study  of  Dynamic  Stability 
Characteristics  of  Para chute -Load  Systems",  Mr,  Edward  J, 
Giebutowski  served  as  Project  Officer, 

The  objective  of  the  effort  was  to  produce  a  computerized 
trajectory  simulation  which  would  describe  the  motion  of  a 
single  narachute  and  its  cargo  from  the  time  of  release  from 
the  aircraft  to  the  time  of  impact. 

‘This  report  describes  the  analytical  modeling  and  design 
rationale  for  the  computer  simulation  of  the  trajectory  and 
shows  typical  results  obtained  from  computer  solutions  of 
several  test  cases  involving  a  selected  number  of  parachute 
configurations  as  used  by  the  US  Army  for  airdrop  of  cargo 
and  personnel. 
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ABSTRACT 


A  method  of  total  trajectory  simulation  was  established 
which  is  based  on  the  governing  equations  of  the  various 
phases  of  an  airdrop  or  recovery  system.  In  view  of  these 
equations,  a  computer  program  capable  of  predicting  the 
performance  characteristics  of  a  parachute -load  system  from 
the  instant  of  initiation  to  the  moment  of  landing  was 
established.  Calculations  were  performed  for  a  number  of 
different  aerial  delivery  systems.  The  calculated  results 
fall  well  within  the  broad  ranges  of  expected  performance, 
based  upon  a  familiarity  with  field  test  results.  The  system 
is  ready  to  be  used  for  overall  prediction  of  parachute  per¬ 
formance  characteristics  and  an  intensive  comparison  of  cal¬ 
culated  and  recorded  field  test  results  is  highly  desirable 
for  validation  and  improvement  of  the  technique  of  total 
trajectory  simulation. 

The  report  is  presented  in  two  volumes.  Volume  I  dis¬ 
cusses  and  shows  the  development  of  the  phases  of  the  total 
trajectory  simulation,  and  includes  the  results  of  sample 
calculations  from  release  to  ground  impact  for  several  cases 
of  parachute- load  systems.  Volume  II  presents  the  computer 
program  and  methods  for  obtaining  numerical  results  in  detail. 
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I,  INTRODUCTION 


The  process  of  delivering  a  load  from  an  airplane  in 
flight  to  the  ground- -or  the  recovery  of  a  manned  or  instru¬ 
mented  capsule,  drone  airplane,  or  missile — encompasses  a 
number  of  dynamic  and  aerodynamic  problems.  The  various 
phases  of  such  a  process  are  usually  organized  in  detailed 
groupings  concerning  parachute  deployment,  inflation  and 
descent.  In  each  field  intensive  investigations  have  been 
carried  out  and  generally  the  basic  physical  and  mathematical 
models  have  been  established.  In  many  cases  the  actual  pro¬ 
cesses  are  so  complicated  that  closed,  or  semi-closed,  solu¬ 
tions  can  be  established  only  by  imposing  stringent  simplifi¬ 
cations  and  restrictions.  In  certain  phases  of  the  entire 
operation,  for  example  the  parachute  deployment  and  inflation, 
the  assumption  that  the  parachute- load  system  moves  in  one 
plane  is  quite  acceptable.  However,  during  the  descent 
phase,  many  types  of  parachutes  impose  motions  upon  the 
system  in  three  dimensions,  and  neglecting  this  fact  is  often 
an  unacceptable  simplification. 

The  availability  of  large  high  speed  computers  makes  it 
possible  to  pursue  numerical  solutions  which  allow  the  intro¬ 
duction  of  non-linear  performance  characteristics  and  consid¬ 
erations  of  multiple  degrees  of  freedom.  This  is  particularly 
important  for  the  analysis  of  the  descent  phase  of  the  opera¬ 
tion  where  large  deflections  of  the  system  from  the  ideal 
trajectory  occur,  and  where  three  dimensional  oscillations 
and  coupling  between  planar  oscillations  and  roll  modes  of 
the  parachute -load  system  can  exist. 

One  objective  of  the  following  study  was  then  the  estab¬ 
lishment  of  a  computer  program  to  calculate  all  performance 
characteristics  and  capable  of  three  dimensional  calculations 
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for  the  descent  phase,  using  established  performance  para¬ 
meters  and  functional  relationships  with  a  minimum  of  numeri¬ 
cal  simplifications.  A  further  objective  was  that  the  pro¬ 
gram  allow  the  continuous  performance  calculation  of  a  given 
parachute- load  system  from  initiation  of  parachute  deployment 
to  landing. 

In  view  of  these  primary  objectives,  mechanical  and 
mathematical  models  for  the  various  phases,  whose  validity 
had  previously  been  proven  by  comparison  with  experimental 
results,  were  selected  from  the  literature  and  formed  the 
basis  for  the  computer  program.  A  major  feature  of  the  pro¬ 
gram  is  that  it  has  been  arranged  so  that  entire  calculation 
methods  can  easily  be  replaced  or  updated  when  better  methods 
become  available.  This  can  be  done  simply  by  replacing  major 
subroutines. 

The  basic  approach  to  the  various  phases  can  be 
described  as  follows.  The  motion  of  the  parachute-load 
system  during  the  phases  of  load  extraction,  when  appi  -.able, 
parachute  deployment,  and  inflation  is  considered  to  be  two- 
dimensional.  The  calculation  follows  the  state-of-th®  art 
methods  and  is  coupled  with  the  equations  for  a  two- 
dimensional  trajectory.  In  the  final  descent  phase,  the 
parachute- load  system  is  treated  as  a  rigid  body  and  a  classi¬ 
cal  mechanics  formulation  is  used  to  find  the  equations  of 
motion  with  six  degrees  of  freedom  which  provide  trajectory 
data  in  three  dimensions,  or  with  three  degrees. of  freedom  for 
two-dimensional  solutions. 

Numerical  results  are  obtained  using  finite-difference 
methods  for  the  trajectory  simulation  prior  to  full  inflation, 
and  the  Runge-Kutta  method  is  used  for  numerical  integration 
of  the  equations  of  motion  for  the  final  descent  phase. 

The  report  is  presented  in  two  volumes.  Volume  I  dis¬ 
cusses  and  shows  the  development  of  the  phases  of  the  total 
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trajectory  simulation,  and  includes  the  results  of  sample 
calculations  from  release  to  ground  impact  for  several  cases 
of  parachute- load  systems.  Volume  II  presents  the  computer 
program  and  methods  for  obtaining  numerical  results  in 
detail. 

The  authors  wish  to  express  their  gratitude  to 
Mr.  Edward  J.  Giebutowski  of  the  United  States  Army  Natick 
Laboratories  for  his  cooperation  in  establishing  the  composi¬ 
tion  and  physical  data  of  the  various  airdrop  systems  and 
for  his  many  constructive  contributions  throughout  this  study. 


II.  SEPARATION  AND  DEPLOYMENT  SYSTEMS 


This  section  describes  the  various  physical  airdrop 
systems  to  which  the  total  trajectory  simulation  applies. 

The  basic  phases  of  functioning  for  any  parachute- load 
system  are:  1)  separation  from  the  aircraft,  2)  deployment 
of  the  parachute,  3)  inflation,  and  4)  free  descent.  The 
physics  of  the  last  two  processes  is  essentially  the  same  for 
any  parachute  type  decelerator;  however,  there  are  many  dis¬ 
tinctly  different  physical  processes  used  to  accomplish 
separation  and  deployment.  This  study  considers  four  different 
systems  which  represent  methods  most  likely  to  be  used  by  the 
United  States  Army  for  cargo  and  personnel  airdrop  (Ref  1). 

The  foui  separation  and  deployment  systems  incorporated 
in  the  trajectory  simulation  are:  1)  static  line,  2)  static 
line  deployed  pilot  chute,  3)  extraction  parachute,  and  4) 
reefed  main  parachute  extraction.  Separation  from  the  air¬ 
craft  is  achieved  by  the  load  merely  falling  away  from  the 
aircraft  in  the  first  two  systems,  whereas  in  the  last  two 
systems  a  parachute  is  attached  to  the  load  to  increase  the 
relative  velocity  between  aircraft  and  load.  A  detailed 
description  of  the  four  systems  follows. 


A.  Static  Line  System 

Figure  1  shows  the  static  line  separation-deployment 
system  used  with  calculations  for  the  T-10  personnel  parachute 
and  the  G-13  cargo  parachute.  The  static  line  is  attached 
to  the  aircraft  while  its  free  end  pulls  the  deployment  bag 
from  the  parachute  after  the  lines  are  extended.  The  canopy 
is  then  free  to  inflate.  No  distinction  is  made  between 
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Static  Line  System 


’’lines  first”  and  ’’canopy  first”  deployments,  since  the 
modeling  equations  apply  to  either. 


B.  Static  Line  Deployed  Pilot  Chute  System 

Figure  2  illustrates  the  use  of  a  static  line  with  a 
pilot  chute,  used  for  separation -deployment  of  the  G-12D 
cargo  parachute.  After  release  of  the  parachute -load  system, 
the  sta:ic  line  deploys  the  pilot  chute  from  its  pack.  The 
pilot  chute  is  then  used  to  deploy  the  main  canopy  in  a 
lines-first  fashion,  and  the  deployment  bag  is  separated. 


C.  Extraction  Parachute  System 

Figure  3  illustrates  the  extraction  parachute  system, 
used  for  separation-deployment  of  the  G-12D  and  G-11A  cargo 
parachutes.  The  load  is  separated  from  the  aircraft  by  means 
of  the  extraction  parachute.  As  the  load  leaves  the  aircraft, 
the  extraction  parachute  force  is  transferred  to  the  deploy¬ 
ment  bag  of  the  main  parachute,  and  the  canopy  is  deployed  in 
a  lines-first  manner.  Figure  4  shows  a  reefing  sequence,  as 
used  with  large  parachutes.  The  extraction  parachute  again 
separates  the  main  canopy  deployment  bag  from  the  load, 
paying  out  the  risers  and  suspension  lines.  But,  as  the 
suspension  lines  tighten,  the  reefing  cutters  are  armed,  and 
the  remainder  of  the  canopy  is  deployed.  After  a  prescribed 
delay  time,  the  reef'.ng  line  is  cut,  and  the  parachute  is  free 
to  inflate  fully  or  to  a  second  reefed  stage. 


o 


Fig  2  Static  Line  Deployed  Pilot  Chute  System 


Extraction  Parachute  System 


D.  Reefed  Main  Parachute  Extraction  System 

The  fourth  system.  Fig  5,  involves  the  use  of  the  main 
canopy  in  a  reefed  configuration  as  the  extraction  parachute. 
The  load  is  pulled  from  the  aircraft,  and  after  a  prescribed 
delay  the  reefing  line  is  severed. 
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ig  5  Reefed  Main  Parachute  Extraction  System 


Ill*  TRAJECTORY  SIMULATION  METHODS 


The  total  trajectory  simulation  was  divided  into  phases 
which  represent  the  significant  physical  processes  of  the 
airdrop  and  relate  to  computer  program  major  subroutines, 
which  can  be  replaced  individually.  Details  of  the  differ¬ 
ent  phases  are  described  in  the  following.  The  first  two 
Subsections,  entitled  Atmospheric  Density  Function  (III.  A) 
and  Two  Dimensional  Point  Mass  Trajectory  Analysis  (III  B), 
represent  cal  relations  which  are  made  at  various  points 
throughout  the  total  trajectory  simulation  and  thus  are 
discussed  first. 


A.  Atmospheric  Density  Function 

During  the  trajectory  simulation,,  the  values  for  the 
air  density  used  in  the  various  equations  governing  the 
motion  throughout  the  airdrop  are  determined  by  the  position 
of  the  parachute- load  system  center  of  mass.  The  function 
that  is  used  to  describe  the  air  density  for  a  standard  day  is 

f>  =  f>o  e'V’/3Z,‘U(’  (  04.  Vi  i  I  S', OOO 

p-  p  (j. 07132.')  /s-,000  ft 

I  r  0  V  ’  <  35,000  ft 

where  p0  «  0.002378  slug/ ft"*.  This  function  was  determined 
by  fitting  two  straight  lines  to  a  logarithmic  plot  of  o 
vs  altitude  determined  from  Ref  2  as  shown  in  Fig  6. 


B.  Two-Dimensional  Point  Mass  Trajectory  Analysis 

The  motion  of  the  parachute- load  system  is  assumed  to 
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Air  Density  Funtion 


be  two- dimensional  during  the  processes  of  separation,  deploy¬ 
ment,  and  inflation,  since  all  forces  act  in  one  vertical 
plane  during  these  processes.  For  the  periods  during  this 
two-dimensional  portion  of  the  airdrop  when  the  position 
of  the  mass  center  does  not  change  significantly,  such  as  in 
separation  and  deployment,  the  motion  of  the  mass  center  is 
approximated  by  a  point  mass  ballistic  trajectory.  Equations 
for  the  trajectory  of  the  parachute -mass  system  can  then  be 
found  by  considering  the  forces  acting  on  a  point  mass  with 
a  given  drag  area. 

Summing  forces  parallel  and  perpendicular  to  the  tra¬ 
jectory  of  a  particle  of  mass  m  with  drag  area  CDS  (Fig  7) 
yields  the  equations 


w>  ^  *  vv^  cosO-lpvtCoS 


(2) 


ae 

at 


W  V  is 

W»v  — —  =  -  w\a  0 


(3) 


The  velocity  and  angle  can  then  be  found  numerically  by 
expressing  Eqns  (2)  and  (3)  in  finite  difference  form,  and 
the  values  of  x  and  z  are  easily  found  from  w  v  sin  0 
and  “  v  cos  ®  expressed  in  finite  difference  form. 
Summarizing,  for  all  processes  which  can  be  approximated  by 
point  mass  trajectories,  the  governing  equations  are 
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C.  Separation  from  Aircraft 

There  are  several  possible  calculation  schemes  for  this 
phase  of  performance,  depending  on  the  separation- deployment 
system  to  be  used.  This  is  the  first  phase  of  the  total 
trajectory  simulation,  beginning  at  the  instant  the  load  is 
released  from  the  aircraft. 


1.  Load  Extraction 

If  the  load  is  to  be  extracted  with  extraction  parachutes 
or  by  the  reefed  main  parachute,  the  governing  equation  is 
for  horizontal  deceleration  unti?  the  load  leaves  the  aircraft, 
i.e. 

-  -  4  Pv*  Ct>S  (8) 


The  modeling  equations  incorporated  in  the  total  trajectory 
simulation  are  then,  expressing  Eqn  (7)  in  finite  difference 
form, 


f  v*C  oS 
2  w> 


(9) 


A  *  »  v  At 


(10) 


where  m  is  the  mass  of  the  entire  recovery  system  and  C^S  is 
the  drag  area  of  the  extraction  parachute,  (C£)Sex)  >  or  the 
reefed  main  parachute,  (On  Sp) .  Equations  (9)  and  (10) 
govern  the  motion  until  the  load  has  traveled  a  specified 
distance  L  and  leaves  the  aircraft,  i.e.,  until 


v0-t  -  *  a  L 
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where  vQ  is  the  aircraft  velocity.  Once  the  conditions  are 
such  that  (11)  is  satisfied,  the  system  is  then  governed 
by  the  two-dimensional  trajectory  Eqns  (4)  through  (7),  with 
the  mass  remaining  the  same  but  increasing  the  drag  area  to 
include  the  load  and  main  parachute  bag  which  are  now  exposed 
to  the  flow,  for  a  specified  time  until  the  extraction  main 
canopy  is  disreefed  or  deployment  of  the  main  canopy  by  the 
extraction  parachutes  is  initiated.  The  next  phase  of  the 
simulation  is  then  inflation  of  the  main  parachute  from  its 
reefed  configuration  for  the  former  case  or  deployment  and 
snatch  force  for  the  latter. 


2.  Static  Line  Extension 

When  a  static  line  is  used,  the  system  is  governed  by 
Eqns  (4)  through  (7)  using  the  mass  and  drag  area  of  the  load 
and  the  packed  recovery  system.  This  is  true  until  the  dis¬ 
tance  between  the  aircraft  and  the  load  equals  the  length  of 
the  static  line,  plus  the  suspension  system  and  half  the 
diameter  of  the  pilot  parf-.hute,  if  used.  This  criterion  is 

VW-xf  +  =  L  Si* lie  +  (12) 

where  L*  is  the  additional  length  of  the  suspension  system 
and  deployed  pilot  parachute  canopy.  If  a  pilot  chute  is 
used,  the  trajectory  can  be  governed  by  Eqns  (4)  through  (7) 
with  the  drag  area  increased  to  include  the  pilot  chute  until 
a  specified  time  is  reached. 


D.  Deployment  and  Snatch  Force 

Once  the  system  has  separated  from  the  aircraft  and 
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coasted  on  the  extraction  or  pilot  parachutes  if  desired, 
i.e„,  remained  unchanged  in  configuration  while  following 
its  trajectory  for  a  specified  time,  the  deployment  of  the 
main  parachute  begins.  Two  methods  are  considered,  static 
line  deployment  and  deployment  by  pilot  or  extraction  para¬ 
chute.  It  has  been  assumed  that  the  "snatch  force"  from  a 
static  line  deployment  can  be  neglected. 


1.  Deployment  by  Static  Line 

For  the  case  of  deployment  by  static  line,  the  trajec¬ 
tory  is  governed  by  the  same  equations  as  during  the  static 
line  separation  process.  Thus  the  governing  equations  are 
Eqns  (4)  through  (7),  and  the  canopy  is  deployed  when  Eqn  (12) 
is  satisfied,  where  L'  is  the  length  of  the  suspension  system 
plus  half  the  diameter  of  the  main  canopy.  The  total  tra¬ 
jectory  simulation  then  proceeds  to  the  inflation  calculation 
phase . 


2.  Deployment  by  Pilot  or  Extraction  Parachute 

The  physical  process  involved  in  this  method  begins  at 
the  instant  that  deployment  of  the  main  parachute  is  initiated 
and  includes  the  deployment  of  the  suspension  lines,  risers, 
etc.,  snatch  force,  and  the  deployment  of  the  main  parachute 
canopy. 

Certain  assumptions  concerning  the  motion  of  the  primary 
and  secondary  bodies  during  this  phase  of  operation  must  be 
made.  Basically,  it  is  assumed  that  the  primary  body  alone 
determines  the  trajectory  that  the  system  follows  during  the 
deployment,  i.e.,  that  the  velocity  vector  of  the  secondary 
body  is  always  aligned  with  that  of  the  primary  body 
(Fig  0) . 


a.  Suspension  Line  and  Riser  Deployment 

The  motion  during  deployment  is  given  by  a  point  mass 
trajectory  analysis  of  the  primary  body  with  the  secondary 
body  decelerating  at  a  different  rate  but  in  the  same  direc¬ 
tion  as  the  primary  body.  To  simplify  the  problem  of  time 
varying  mass  as  the  suspension  system  is  deployed,  the  pri¬ 
mary  and  secondary  body  masses  are  defined  in  a  manner  simi¬ 
lar  to  Ref  3  during  deployment  of  the  suspension  system, 

v*  -  Wj  *  \  "V,  (13) 

Mx  =  +  4^33  4  '"'fb  (14) 


where  m_  is  the  mass  of  the  suspension  lines,  risers, 
s  s 

extensions,  bridle,  and  links  and  m^  is  the  mass  of  the  pilot 
or  extraction  parachute  and  the  main  parachute  bag.  In  finite 
difference  form,  the  governing  equations  during  deployment 
of  the  suspension  system  are  then 


A0  =  -  At 


(15) 
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Avx=  (3co,0-  )  At  (16) 

Av1I=  C3t.se-  )  At  (17) 

4*  =  vz  cos  6  At  (18) 
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Vx  COS  0  &i 


*1  =  V*At-V*A t  (20) 

where  i  is  the  distance  between  the  primary  and  secondary 
bodies.  C^Sjj  in  Eqn  (17)  is  the  drag  area  of  the  pilot  or 
extraction  parachute  and  the  main  parachute  bag. 

b.  Snatch  Force 

When  the  value  of  X ,  which  is  initially  zero,  equals 
the  length  of  the  suspension  lines,  risers,  extensions,  and 
bridle,  the  line  elongation  and  the  development  of  the  line 
force  begins.  The  force  calculation  follows  the  energy 
conservation  method  of  Ref  3,  and  the  maximum  or  snatch  force 
is  given  by 
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V  =  MIV1  *  (27) 

and  A  is  the  inverse  of  the  effective  spring  constant  of  the 
suspension  lines,  risers,  extensions,  and  bridle.  The  value 
Mj  is  given  by  +  ngs,  and  Eqns  (26)  and  (27)  reflect 

the  assumption  that  the  suspension  system  is  moving  at  the 
same  velocity  as  the  load  and  only  the  main  parachute  canopy 
is  accelerated  at  snatch.  The  time  increment  during  which 
the  lines  elongate  is  very  short  and  is  neglected  in  view  of 
the  total  elapsed  time. 

c.  Main  Parachute  Canopy  Unfolding 

After  snatch,  the  main  parachute  unfurls  from  the  bag. 

The  primary  body  is  now  considered  to  consist  of  the  suspended 
load,  risers,  and  canopy,  and  the  secondary  body  consists  of 
the  pilot  or  extraction  parachute  and  the  deployment  bag. 

The  equations  governing  the  motion  during  unfolding  of  the 
main  parachute  are  similar  to  Eqns  (14)  through  (19)  with  the 
exception  that  (15)  and  (16)  are  replaced  by 
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The  drag  area  in  Eqn  (28)  is  increased,  by  the  addition  of 
1.5%  of  the  fully  inflated  drag  area  of  the  main  parachute. 
This  value  represents  the  drag  area  of  the  uninflated  main 
parachute  and  was  taken  from  Ref  4.  The  trajectory  is  cal¬ 
culated  with  Eqns  (15),  (28),  (29),  (18)  and  (19)  and  the 
separation  with  Eqn  (20)  until  the  value  of  ji  equals  the 
value  at  snatch  plus  half  the  nominal  diameter  of  the  main 
parachute  at  which  time  the  canopy  is  deployed.  The  next 
phase  is  then  the  inflation  phase. 


3.  Summary  of  Assumptions  for  the  Deployment  Phase 

In  the  following  the  assumptions  on  which  the  preceding 
calculations  were  based  are  summarized.  For  static  line 
deployment,  it  is  assumed  that  (1)  the  static  line  remains 
in  a  straight  configuration  between  the  aircraft  and  the 
load,  and  (2)  there  is  no  appreciable  "snatch  force."  When 
the  main  decelerator  is  deployed  by  extraction  or  pilot  para¬ 
chute  several  assumptions  are  made:  (1)  the  trajectory  of 
the  parachute -load  system  during  deployment  is  determined 
only*by  the  load;  (2)  the  positions  of  the  primary  and  second¬ 
ary  bodies  can  be  described  by  point-mass  trajectory  equa¬ 
tions  using  "average"  values  for  the  masses;  (3)  the  occur¬ 
rence  of  snatch  is  instantaneous  for  elapsed  time  consideration; 
(4).  the  mass  of  the: .‘parachute  bag  and  pilot  or  extraction 
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parachute  is  decelerated  relative  to  the  load  as  the  canopy 
unfurls.  The  above  assumptions  are  simplifications  of  actual 
deployment  conditions.  The  need  for  a  more  exact  treatment 
of  this  phase  should  be  determined  in  future  work. 


E.  Inflation  of  the  Main  Parachute 

The  inflation  of  the  main  parachute  begins  at  the  instant 
the  parachute  is  deployed  in  a  stretched-out  configuration 
with  assumed  zero  projected  area,  or  at  some  initiation  time 
tjj  when  the  parachute  is  used  in  a  reefed  configuration  to 
extract  the  load.  The  total  trajectory  simulation  calculates 
the  motion  of  the  parachute -load  system  as  it  inflates  to  its 
fully  open  configuration  with  or  without  intermediate  reefed 
configurations. 


1.  Governing  Equations 

The  basis  for  the  total  trajectory  simulation  during  the 
inflation  of  the  main  parachute  is  a  coupling  of  the  calcula¬ 
tions  of  the  opening  force  (Ref  5),  and  a  ballistic  trajectory. 
The  opening  force  calculation  assumes  that  the  inflating 
parachute  has  the  shape  of  a  truncated  cone  with  a  hemispheri¬ 
cal  cap  such  that  the  suspension  lines  extend  parallel  to  the 
edges  of  the  truncated  cone  (Fig  9).  The  trajectory  calcula¬ 
tion  is  based  on  force  equations  parallel  and  perpendicular 
to  the  trajectory  of  the  mass  center  of  the  parachute- load 
system,  and  it  is  assumed  that  the  velocity  components  of 
the  load  and  the  parachute  normal  to  the  systems  axis  are 
assumed  to  be  negligible.  This  assumption  may  be  questionable, 
especially  for  large  parachutes  with  long  distances  between 
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Parachute -Load  System  During  Inflation  of  the 
Main  Parachute 


parachute  and  load;  however,  it  was  felt  that  such  an 
assumption  was  a  reasonable  first  approximation.  It  should 
be  noted  that  the  above  velocity  components  are  not  neglected 
in  the  free  descent  phase. 

Therefore,  the  parachute -load  system  is  considered  to 
be,  in  effect,  a  point  mass  with  varying  mass  and  drag  area 
during  the  inflation  process.  The  force  equations  parallel 
and  perpendicular  to  the  trajectory  of  the  mass  center  are 
then 
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the  equations  governing  the  motion  can  be  written  in  finite 
difference  form  as 
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Thus  the  motion  of  the  parachute -load  system  during  infla¬ 
tion  can  be  found  numerically  from  Eqns  (33)  and  (34)  if  the 
values  of  m . ,  ma,  dm. /dt,  dm  /dt,  and  0  are  known  throughout 

1 9  fl 9  1  u  P 

the  inflation  process. 

a0  inflation  Without  Reefing 

The  basic  assumptions  of  Ref  5  are  that  the  canopy 
geometry  at  any  time  instant  during  inflation  is  determined 
by  a  linear  growth  of  the  projected  area  and  a  linear  function 
of  the  influx  versus  a  dimensionless  time  T  *  (t  -  tQ)/t 
These  are  commonly  used  assumptions  which  are  subject  to 
change  on  the  basis  of  future  studies.  Certain  differences 
appear  between  the  definitions  of  filling  time  used  in  this 
report  and  in  other  literature,  such  as  Ref  6.  The  dimension¬ 
less  time  scale  T  as  used  in  Ref  6  ranges  from  T  »  0  at  peak 
snatch  force  to  T  =  1  when  the  projected  area  first  reaches 
its  steady  state  value..  The  dimensionless  time  scale  in  this 
study  begins  after  the  main  parachute  canopy  unfolding  process 
has  been  completed,  T  »  0,  and  continues  until  the  steady  state 
volume  of  the  parachute  canopy  is  reached,  T  *  1.  Thus  the 
base  of  the  dimensionless  scale  is  referred  to  as  the  final 
filling  time,  tff. 

The  first  assumption  provides  for  the  projected  diameter 
the  relationship 

Op-  S-T*  <3» 

All  other  terms  in  the  equations  of  motion  can  then  be 
determined  from  Dp  and  the  assumed  geometry.  Expressing 
projected  diameter  as  a  function  of  dimensionless  time  adds 
the  unknown  t£f  to  the  problem.  The  solution  for  t^.  is 
accomplished  by  examining  the  continuity  equation  describing 


the  imbalance  of  mass  flow  into  and  out  of  the  parachute  canopy. 
This  equation  is,  in  accordance  with  Ref  5, 

g  =  7W44V  [ll+2.2cT-T)  £  -  'Jjg*  ]  (36) 

which  is  based  on  the  inflow  velocity  in  the  form  of 

^»(l+2.acT-T)  (37) 

Thus  the  filling  time  tff  for  inflation  without  reefing  can  be 
found  by  integration  of  Eqn  (36)  as  T  ranges  from  0  to  1, 

The  volume  V  of  the  fully  inflated  parachute  is  equal  to 

V-  n-fc*  S^lO+ZScT-T)^  -  dJ  (38) 

The  filling  time  is  readily  found  by  a  numerical  pro¬ 
cedure  as  follows.  An  initial  value  must  be  estimated  for 
the  filling  time.  Equation  (38)  can  then  be  numerically 
integrated  with  Simpson's  Rule  using  a  number,  n,  of  equal 
increments  of  time  AT  such  that  T  ranges  from  0  to  1,  yielding 
a  result  which  depends  on  the  assumed  filling  time.  This 
numerical  integration  is  expressed  as 

['V + 2 vt ♦  4  t»a*2v< ♦  •  ••  +  vjj  (39) 

where 


Tj,  j  =  1,  2,  ...,  n,  represent  the  values  of  the  time  scale 
selected  for  the  numerical  integration.  The  numerical  value 
of  D  (T.)  is  found  from  Eqn  (35)  evaluated  at  T  13  T .,  and 

i  J  J 
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the  corresponding  numerical  value  of  d(Tj)  is  determined  by 
the  geometrical  relation 

4U.*LOop 


2  0,-*  Of 


(41) 


The  velocity  v(T.) ,  however,  must  be  found  from  the  equations 
of  motion,  Eqns  (33)  and  (34)  with  initial  conditions  9=  90 
and  v  -  v  for  T  =  0.  The  mass  terms  in  the  equations  of 
motion  based  on  the  assumed  geometry  and  expressed  as  functions 

of  T  and  D  are 

up  Op 

(42) 
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gnp  d.4 

(Pp«»/)*  <*-T 


4  p  ‘ 0f 

I’S  ,  ^p-T-p 
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where 


_  _Do_  (46) 

dT  "  -K1'1* 


Eqns  (32)  and  (33)  are,  in  terms  of  tj£  and  T, 


Thus  the  velocity  can  be  found  for  any  value  of  Tj  by  using 
Eqns  (48)  and  (49)  for  the  particular  choice  of  t^,  and  the 
numerical  integration  (39)  can  be  performed. 

The  value  of  V(tffj)  determined  from  Eqn  (39)  is  then 
compared  with  the  known  volume  of  the  fully  inflated  canopy, 
and  a  new  approximation  to  the  filling  time  is  found  from 
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V 


(50) 


The  process  outlined  above  is  repeated  until  jV  -  V(tff^) 
is  small  enough  that  the  filling  time  can  be  considered 
accurate , 

It  may  be  mentioned  that  the  method  of  Ref  5  follows  a 
similar  scheme,  however,  the  various  functions  are  simplified 
by  means  of  curve  fitting  since  a  practically  closed  solution 
is  pursued  in  Ref  5* 

Once  the  filling  time  is  known,  the  trajectory  of  the 
parachute- load  system  mass  center  is  determined  from  Eqns 
(48)  and  (49)  in  conjunction  with 


Ax  a 

v  -fc^  AT 

(51) 

A*- 

V  cos  0  atyr  AT 

(52) 

and  the  instantaneous  opening  force  is  found  from 


b,  Reefed  Inflation 

When  the  inflation  process  includes  intermediate  reefed 
configurations,  the  trajectory  simulation  is  divided  into 
periods  of  inflation  and  periods  of  coasting  in  the  reefed 
configuration,  if  desired.  Provision  for  the  coasting 
periods  is  made  by  prescribing  a  reefing  cutter  delay  from 
the  time  of  suspension  line  deployment,  or  from  t  *  0  when 
the  load  is  released  in  the  aircraft  for  systems  extracted 
by  the  reefed  main  parachute. 
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During  the  coasting  periods,  the  drag  area  remains  con¬ 
stant,  and  the  motion  is  described  by  the  two-dimensional 
ballistic  trajectory  equations  (4)  through  (7),  with  the  value 
of  the  load  drag  area  plus  the  drag  area  corresponding  to  the 
particular  reefed  configuration  of  the  parachute,  and  the 
mass  of  the  parachute- load  system.  When  the  final  configura¬ 
tion  of  an  inflation  phase  is  the  fully  inflated  configuration, 
there  is  no  coasting  period,  but  the  trajectory  simulation 
proceeds  immediately  with  the  free  descent  phase. 

During  each  period  of  inflation,  it  is  assumed  that  the 
projected  area  increases  linearly  with  time  (Fig  10) , 

Defining  a  dimensionless  time  scale  for  each  period  of 
inflation  such  that 


(54) 


where  t  is  the  initial  time  and  tf-  is  the  inflation  time 
o 

for  the  particular  reefed  inflation  period,  enables  the 
problem  to  be  viewed  in  a  manner  similar  to  the  approach 
followed  for  the  unreefed  inflation.  The  projected  area  growth 
for  each  of  the  inflation  periods  is  given  by  (Fig  10) 


(55) 


where,  with  the  geometry  as  used  before, 

Ujl  <KLs-‘L(0R(>+ 
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(56) 
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In  these  equations,  R 0  and  are  the  reefing  ratios  defined 
by 


(57) 


corresponding  to  the  reefing  lines  at  the  beginning  and  end 
of  the  inflation  period.  The  projected  diameter  is  then 
found  from  Eqn  (55)  to  be 

Dp=  O.Vv'.*Tr  ♦  C I-TrV  (58) 


Comparing  Eqns  (35)  and  (58),  it  is  evident  that  one  can 
consider  the  projected  diameter  growth  during  an  inflation 
period  to  be  a  function  of  T,  the  dimensionless  time  scale 
which  would  exist  if  the  inflation  were  taking  place  in  one 
complete  step  without  reefing,  i.e.,  the  functions  describing 
the  reefed  inflation  can  be  considered  to  be  portions  of  the 
functions  developed  for  the  unreefed  case.  The  values  of 
projected  diameter,  masses,  etc.  can  then  be  found  for  any 
value  of  TR  by  evaluating  the  equations  presented  in  Item  a 
at  a  value  of  T  given  by 

T *  if  +  Vi*  (J-Tr)'  (59> 


which  is  obtained  by  equating  the  right  hand  sides  of  Eqns 
(35)  and  (58). 

The  procedure  for  calculation  of  the  filling  time,  tra¬ 
jectory,  and  opening  force  is  similar  to  that  discussed  for 
the  unreefed  case.  The  filling  time  is  found  from  numerical 
integration,  with  TR  ranging  from  0  to  1,  of 

V* =  -K-fc+ft  v  l  0*  M  cT-T)  £  -  LifcSe  ]  (60) 
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i.e 


+  (61) 

where  Vj  is  given  by  Eqn  (40)  with  Tj  evaluated  from  Eqn  (59). 
The  value  represents  the  increase  in  the  volume  of  the 
parachute  canopy  *s  it  inflates  from  the  initial  to  the  final 
reefed  configuration  and  is  given  by 


Successive  values  of  tf  are  given  by 

*i 


until  the  value  J V(t^R  )  -  V^j  is  small  enough. 

The  necessary  values  for  evaluation  of  Eqn  (61)  are 
determined  as  follows. 


The  velocity  is  found  from 


*  v*  - 
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The  values  of  Dp,  m^,  ma,  dm^/dT,  and  dma/dT  are  found  from 
Eqns  (35) ,  (42) ,  (43) ,  (44) ,  and  (45)  with  T  from  Eqn  (59) . 
The  increment  AT  is  determined  from  the  prescribed  increment 

iTR  by 

AT=  AT  (66) 

The  trajectory  is  described  by 

Ay  =  v  s\v->  $  ATg  (67) 

As*  v  cos  0  -fc^  ATr  (68) 


and  the  instantaneous  opening  force  is  given  by 

<V  ■£-«*)  l69> 


2.  Summary  Assumptions  for  the  Inflation  Phase 

The  inflation  phase  of  the  trajectory  simulation  is 
modeled  by  the  following  assumptions.  (1)  The  shape  of 
the  parachute  is  described  by  a  truncated  cone  with  a  hemi¬ 
spherical  cap  and  the  projected  area  increases  linearly  with 
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time  during  both  reefed  and  unreefed  inflation  periods. 

(2)  The  motion  is  confined  to  a  vertical  plane.  (3)  The 
magnitudes  of  the  components  of  parachute  and  load  veloci¬ 
ties  normal  to  the  trajectory  of  the  mass  center  are  assumed 
negligible.  (4)  When  reefing  is  used  and  coasting  between 
the  inflation  periods  is  arranged,  the  path  of  the  center  of 
mass  can  be  approximated  by  the  two-dimensional  point  mass 
trajectory  equations.  (5)  The  time  functions  describing 
the  performance  of  the  parachute  when  inflating  without  reef¬ 
ing  can  be  used  to  describe  the  reefed  inflation  periods  by 
adjusting  the  time  scale. 


F.  Free  Descent 

This  phase  begins  at  the  instant  the  parachute  is  fully 
inflated  and  terminates  at  ground  impact. 

Two-body  trajectory  analyses  are  available,  such  as 
Refs  7  and  8,  which  describe  the  two-dimensional  motion  of 
the  separate  masses  of  the  parachute  and  payload  when  con¬ 
strained  by  the  elastic  connection  (risers,  slings,  and 
suspension  lines)  between  them.  While  consideration  of 
elasticity  is  undoubtedly  a  more  accurate  representation  of 
an  actual  system,  it  was  felt  that  a  reasonable  approximation 
of  the  free  descent  phase  could  be  achieved  by  analyzing  the 
system  as  a  rigid  body  having  realistic  distribution  of  the 
masses.  It  should  be  noted  that  the  systems  under  consideration 
are  limited  to  single  parachutes  with  assumed  point  mass 
payloads.  Further,  the  requirement  for  developing  three- 
dimensional  equations  of  motion  and  computer  programs  would 
have  been  considerably  complicated  by  the  assumption  of  a 
two- body  system. 
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1.  Governing  Equations 

In  the  following,  a  set  of  orthogonal  coordinate  axes 
fixed  in  the  parachute  load  system  is  defined  (Fig  11a)  with 
respect  to  which  the  mass  distribution  is  constant.  These 
body- fixed  axes  are  denoted  by  XYZ.  The  Newtonian  or  space - 
fixed  axes,  with  respect  to  which  all  accelerations  must  be 
measured,  are  chosen  to  be  fixed  at  the  point  of  release  from 
the  drop  aircraft.  These  Newtonian  coordinates  are  the  same 
axes,  denoted  xyz,  used  to  describe  the  motion  of  the  system 
during  the  previous  phases  of  the  trajectory  simulation. 

The  orientation  of  the  parachute -load  system  is  described 
by  Euler  angles,  f,  0,  cp,  chosen  to  correspond  with  the 
notation  of  Refs  9  and  10.  The  Euler  angles  chosen  for  the 
trajectory  simulation  are  three  consecutive  rotations,  begin¬ 
ning  with  a  coordinate  system  x^y^z^  parallel  to  the  space- 
fixed  reference  frame  (see  Fig  lib),  defined  as  follows: 

(1)  a  rotation  if  about  axis  0Z^,  yielding  ^2^2Z2 

(2)  a  rotation  0  about  axis  0y2,  yielding  x^y^z^ 

(3)  a  rotation  {p  about  axis  0X^,  yielding  XYZ. 

Vector  quantities  expressed  in  the  space -fixed  axes  are  then 
related  to  quantities  expressed  in  the  body-fixed  axes  by  the 
direction  cosine  matrix,  i.e. 


e  -  AE  (70) 

where  e  and  E  are  vectors  in  the  space-fixed  and  body-fixed 
coordinate  systems,  respectively.  The  direction  cosine  matrix 
is  the  direction  cosine  matrix  of  Ref  10  transposed: 
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Fig  11b  Euler  Angles 
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Choosing  the  origin  of  the  body-fixed  coordinate  system 
as  the  mass  center  of  the  parachute-load  system  allows  the 
equations  of  motion  in  the  body  fixed  system  to  be  written  as 
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At- 


(72) 
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(73) 


Time  derivatives  of  vectors  in  the  rotating  reference  frame 
must  be  measured  relative  to  the  Newtonian  reference  frame. 
If  the  parachute,  and  thus  the  body-fixed  coordinate  system, 
rotates  with  an  angular  velocity  jaj,  the  time  derivative  of  a 
vector  E  is  given  by  Euler* s  equation, 

•  ^n,***!  <»> 

Thus  the  equations  of  motion  (73)  and  (74)  can  be  expressed 
in  scalar  form  as  follows: 

t 

f. 
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(75) 


(76) 


Fy  =  +vnp'^['V+  t?\J-  PW] 


£PZ-  [wi  +  w,,J+voj)['W+  PV-(i\j] 


(77) 


z  l=  ~  QIxv  “  ■*■  ^R(^»B“Xvr) 

+C^-©l')Xys  +  PR^xr  -  PC)Txa 
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Gilyy-Rlr,  -  PI xv  +  PRI^*e‘I«) 


+CP,-Rt)lsa  +  PQ^y*  ~  ^Zvc 


(79) 


»  f?lM-  P^xi  -QIy*  +  PQ^Iyy-Ixx)  (g0) 
+  {?*-  Ql)  IXY  +  QR  I*,  -  PRI*, 


where  L,  M,  N  are  the  body-fixed  components  of  the  moment  M, 
U,  V,  W  are  the  body- fixed  components  of  the  center  of  mass 
velocity  v,  and  P,  Q,  R  are  the  body-fixed  components  of  the 
angular  velocity  j.  In  terms  of  the  angular  velocity  com¬ 
ponents,  the  Euler  angles  are  found  from  the  differential 
equations  (Ref  8) 
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Cp  —  P*  Qs^cp-Uv^Q  +  Row  Cp-tc^Q 


(82) 


t  *  ^Q$wi(p  +  RcoiCy')  St c.  0 


(83) 


The  velocity  of  the  mass  center  in  space- fixed  coordinates  is 
related  to  the  velocity  in  body-fixed  coordinates  by 


XsUau  +  Ya,t  +  Woi*3 

(84) 

^  =•  "0  <*21  +  "V < Xtz  ¥  VJ (*23 

(85) 

(86) 

Thus  Eqns  (75)  through  (86)  represent  a  complete 
description  of  the  system,  and  a  solution  to  these  twelve 
equations  provides  the  necessary  information  for  a  three- 
dimensional  trajectory  simulation.  Linearizing  the  equations 
as  is  usually  done  in  dynamic  stability  analyses,  for 
example  airplane  dynamic  stability,  is  not  an  appropriate 
solution  method  for  a  trajectory  simulation  of  a  parachute - 
load  system  due  to  the  large  angle  deflections  which  occur  in 
parachute- load  system  trajectories.  Thus  the  approach  taken 
for  this  trajectory  simulation  is  to  solve  the  equations 


numerically  using  the  Runge-Kutta  method.  The  remaining 
task  which  must  first  be  accomplished  is  the  formulation  of 
the  external  forces  and  moments  in  the  equations  of  motion 
(75)  through  (80) . 

a.  Three  Degrees  of  Freedom 

When  the  trajectory  simulation  is  restricted  to  three 
degrees  of  freedom,  the  quantities  V,  P,  R,  cp,  ijr,  and  y  are 
all  constrained  to  zero.  Figure  12  illustrates  the  orienta¬ 
tion  of  the  parachute -load  system  and  the  external  forces 
and  moments  which  act  upon  the  parachute  and  the  load  for  the 
three-degree-of-freedom-problem.  The  formulation  of  the 
equations  of  motion  is  simplified  considerably  by  the 
restriction  to  three  degrees  of  freedom,  and  phe  number  of 
equations  required  for  a  complete  description  of  the  system 
is  reduced  from  twelve  to  six.  These  equations  are: 


(87) 

(88) 

£  hA  *  Q.  IfY 

(89) 

©  =  Q 

(90) 

X  *  co s  0  Vj  $5v«i  0 

(91) 

2  a  -X7  jivi  0  +  Wcw  0 

(92) 
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The  external  forces  and  moments  which  are  considered  in 
the  trajectory  simulation  are  due  to  the  weights  of  the 
various  components  of  the  parachute- load  system,  the  aero¬ 
dynamic  tangent  and  normal  forces  and  aerodynamic  moment  due 
to  the  parachute,  the  aerodynamic  drag  force  on  the  load,  and 
the  force  due  to  the  acceleration  of  the  included  mass  and 
the  effect  of  the  apparent  mass  of  the  parachute.  The  lengths 
if,  L 23  and  JL ^  from  the  system  mass  center  to  the  load, 
parachute  center  of  volume,  and  moment  center,  respectively, 
are  signed  quantities,  the  positive  direction  coinciding  with 
the  positive  Z-coordinate  direction.  These  distances  are 
functions  of  the  air  density,  since  the  mass  center  position 
changes  as  the  included  mass  changes.  The  aerodynamic  normal 
force  and  the  aerodynamic  moment  due  to  the  parachute  act 
at  the  moment  center,  which  is  defined  as  a  point  located  a 
distance  equal  to  Dq  ahead  of  the  canopy  skirt.  Thus  the 
lengthjt^  is  always  positive,  JL ^  is  always  negative,  and  JL^ 
can  conceivably  have  either  sign. 

The  force  required  to  accelerate  the  included  mass  and 
the  effect  of  the  apparent  mass  are  given  by 

F=  (93) 

assuming  that  the  apparent  mass  has  the  same  value  in  each 
of  the  coordinate  directions.  The  acceleration  of  the  para¬ 
chute  center  of  volume  is  given  by 

(.ISt'lxYZ  +  (94) 

The  components  of  the  forces  for  the  two-dimensional  case 
can  then  be  written  as 
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(95) 


2  Fx  *  ♦  Fw  +  Dx  3mo<x 

-K^^iu+isw)  +  GUx] 

£FZ  »  ^VAjj+^+v^p)  ^c»C  -T  -  Dx  cosolx 
“ (^1+  ^ «.') [iw -  ) -  Q2ltJ 

If  the  moment  of  inertia  tensor  is  defined  such  that  it 
includes  the  effect  of  the  apparent  and  included  mass,  i.e. 

I  -  (97) 


\ 


the  moment  due  to  the  apparent  force  need  not  be  considered. 
The  moment  in  the  two-dimensional  case  is  then 

£M=  *  DX  t,  ♦  (98) 


The  masses  of  the  risers,  suspension  lines,  bridle,  etc.,  are 
neglected  in  the  formulation  of  the  moment  (98) . 

Using  the  definition  of  total  mass  employed  in  the 
inflation  analysis 

yyVm  **1+  Vv>$5«-  vnp  +  4  (99) 


the  equations  of  motion  (87)  through  (89)  can  be  written 
in  terms  of  equations  for  the  derivatives  of  the  velocity 
and  angular  velocity  components  in  the  body-fixed  coordinate 
system.  This  formulation  is 
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(101) 
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The  six  equations  (100),  (101),  (102),  (90),  (91),  and  (92) 
then  form  a  system  of  first  order  differential  equations  which 
completely  describe  the  parachute -load  system. 

The  aerodynamic  forces  and  moments  are  given  by 


T "ip  Gre  So 

(103) 

FN='4p"'ft  Ck i0S0 

(104) 

\  pvp*  D» 

(105) 

D*  =  c0SA 

(106) 

The  velocities  of  the  load  and  parachute  are  found  from 

Vj^  »  (107) 
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(108) 


Vp*.  V1  +  (O  +  <3  laf 


and  the  angles  of  attack  of  the  parachute  and  load  are 


(109) 


(110) 


The  angle  is  required  since  the  coefficients  Cp£q9  C^,  and 
are  dependent  on  the  parachute  angle  of  attack. 

The  lengths  1 2*  i  3  are  found  from  the  geometry 
of  the  parachute -load  system.  Figure  13  shows  the  profile 
of  a  parachute -load  system  in  the  XZ-plane,  which  is  defined 
as  the  plane  of  symmetry  for  a  general  parachute- load 
system.  The  distances  Sp  S2,  s^,  s^,  and  s^  are  reference 
distances  from  the  canopy  skirt  to  the  cr-.ters  of  mass  of 
the  suspension  lines,  risers,  riser  extensions,  bridle,  and 
load,  respectively,  and  sc  is  the  distance  to  the  parachute 
center  of  volume.  If  all  these  distances,  as  well  as  the 
volume  V  of  the  inflated  parachute  canopy  and  the  masses  of 
the  suspension  lines,  risers,  riser  extensions,  bridle,  and 
load  are  knovm,  the  parachute -load  system  mass  center 
location  can  be  determined  for  any  value  of  the  air  density. 
The  distance  of  the  mass  center  from  the  canopy  skirt  is 
then  given  by 


l*  s’*'*'*?  +>rrvt 
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Where 


Vv^  m  pV 

(112) 

It  is  assumed  that  the  values  of  Sp  S2,  s^, 
do  not  change  during  the  free  descent  phase, 
lengths  from  the  system  mass  center  are 

s/s  sc.  and  s 

4*  5*  c 

Thus  the 

(113) 

it-  "S-St. 

(114) 

l.j=  Oe-  s 

(115) 

An  approximation  to  the  moment  of  inertia  of  the  parachute¬ 
load  system  about  the  Y  axis  is 

(116) 

where 

1y  -  vvij  !L^  +  vnLi  ts,-S  )Z  +  w^(S£-  S 

+  Vr»B  4  +•  if  (117) 

and 

^  l*  (118) 
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The  relation  (118)  is  based  on  Ref  11  with  the  coefficient 
0.252  selected  for  a  solid  flat  circular  parachute.  The 
apparent  mass  is  assumed  to  be  given  by 


W*w  T 


(119) 


which  results  from  Eqns  (43)  and  (44)  when  Dp  “ 

Thus  all  information  is  available  for  the  solution  of 
the  equations  governing  the  motion  of  the  jystem  during  the 
free  descent  phase  of  the  trajectory  simulation.  Once  the 
equations  are  solved  for  U,  W,  Q,  0,  x,  and  z,  further 
information  can  be  readily  obtained.  In  particular,  the 
position,  velocity,  and  acceleration  of  the  load  relative  to 
the  space- fixed  coordinate  system  are  given  by 


riy "  *+  0  St>n  0 


(120) 


'(je,  =  A,co*6 


(121) 


vix  l.^cose  +  w  e;v>9  (122) 

vJ,=-(y+(H  +  WcosO  (123) 

alx  (124) 

C<T-*Q"W+Gl».i')siv>©+  (125) 
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b.  Six  Degrees  of  Freedom 

For  the  general  case,  allowing  six  degrees  of  freedom 
for  the  parachute- load  system  during  the  free  descent  phase, 
the  equations  describing  the  system  are  far  more  complicated 
than  those  for  three  degrees  of  freedom.  The  external 
forces  and  moments  are  indicated  in  Fig  14. 

In  order  to  express  the  forces  and  moments  in  the  body- 
fixed  system  it  is  convenient  to  define  the  following  angles: 


w  -  ^  Qgi 


(126) 


Pi*  ^ 


(127) 


17+ Ql- 


(128) 


fV  **" 


(129) 


j  V-  PJt, 


(130) 


k  55  w' 


(131) 


T7*-QJL, 


[iv-piO4  *V],/2- 


where  the  subscripts  1  and  p  refer  to  the  load  and  the  para¬ 
chute  center  of  volume,  respectively.  The  angles  ctf,  and 
are  the  same  angles  of  attack  in  the  XZ-plane  as  discussed 
previously.  The  angles  0^  and  0p,  analagously,  are  angle* 
of  attack  in  the  YZ-plane.  The  angles  and  6 ^  represent 
the  angles  between  the  load  velocity  and  the  XZ-  and  YZ- 
planes,  respectively. 

The  force  required  to  accelerate  the  included  mass  and 
the  effect  of  the  apparent  mass  are  given  by  Eqn  (93)  and  the 
forces  and  moments  are 


*FX=  t’-’V  ™s«,-vr>p')  3  0-3,  +  FWx 

+  Pt.  COi  -  (vo;,  *■  •  (132) 

-  RV*  XAGt^PR.)] 


*  fy  =  ^ 3*3 z  *  fnv 

“D^COiSji  %hn  vw\qJ^  • 

[v+  RV  -  PW-  A2(P-Q*)] 


(133) 


2  F«  =  (vv> 


ljL+vosi  +  w1p)  ga33  -T  -  Da 

- [ w  +  pv -  Q' V  .  i2(p*«. Q1) ] 


(134) 
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v^t,)  _  pwy 

+  i,  COS  6jlSiv)  (Sj_  + 


(135) 


3  Az  *vv'x1-i)  *  Fwx  A3 

+  A,  tosY^  s'>n +  May 


(136) 


£N  =  M& 


(137) 


where  Mc  is  a  control  or  turning  moment  applied  to  the  system. 
Arranging  the  equations  of  motion  as  equations  for  the  time 
derivatives  of  U,  V,  W,  P,  Q,  and  R,  we  have 


v*  (  *  ^  r2£-)  3**  *  *  %  cojVjl  im*K 

-  "V™*  JlAa*PO  -  OVJ  +  «V 

W>T 


(138) 


(139) 


— ^  Jl^P-QR.)  -  R\7  +  PW 


making  use  of  the  fact  that  all  parachutes  are  bilaterally 
symmetric  and  hence  (choosing  the  XZ-plane  as  the  plane  of 
symmetry)  1^  “  IYZ  “  0.  For  most  parachutes,  it  is  also 
true  that  the  XYZ-axes  are  the  principal  axes,  so  that  Ixz  = 
0,  and  1^  -  Iyy*  Equations  (141)  through  (143)  are  left  in 
the  general  form,  however,  and  are  thus  applicable  to  any 
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parachute.  Equations  (137)  through  (142)  and  (81)  through  (86) 
then  form  a  system  of  twelve  differential  equations  which 
can  be  solve'9  for  U,  V,  W,  P,  Q,  R,  0,  cp,  x,  y,  and  z, 

thus  completely  describing  the  velocity  and  angular  velocities 
as  well  as  the  Euler  angles  and  location  of  the  parachute¬ 
load  system. 

The  parachute  aerodynamic  forces  and  moments  and  the 
drag  of  the  load  are  given  by 

T-4fN/p*cTo  S0 

(144) 

(145) 

FmY*  "ip  Vf 

r  * 

(146) 

• 

(147) 

M*?  *  i  f  v(>2  C’MYo  S, 

(148) 

Djj  =  if  VCoS* 

(149) 

where  ? 

(150) 

Vp1  <U y  4-  ^7 -  Pit?  v  Wl 

(151) 
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Again,  the  coefficients  in  Eqns  (144)  through  (148)  depend  on 
the  angles  of  attack  and  (3^. 

Since  the  XZ- plane  is  a  plane  of  symmetry,  the  lengths 
ip  i 2>  an^  value  of  Iyy  are  found  in  the  same 
manner  as  for  three  degrees  of  freedom,  with  Eqns  (113) 
through  (116).  In  general,  the  values  of  1^,  and  1^  are 

not  given  by  such  a  simple  approximation,  but  depend  on  the 
particular  parachute -load  system.  As  discussed  previously, 
it  is  often  the  case  that  1^  =  Iyy  and  *xz  *  but  even  in 
this  case  the  value  of  I ^  must  be  formulated  for  the  general 
six  degree  of  freedom  system. 

Once  the  equations  have  been  solved  for  U,  V,  W,  P,  Q, 

R,  0,  tpf  x,  y,  and  z,  which  are  related  to  the  mass  center, 
the  position,  velocity,  and  acceleration  of  the  load  are 
given  by 


Vjjx  =  x+  A,  a,j  (152) 

Yiy=  (153) 

(154) 

\7<x,s  (155) 

i)  Ottt  +  tv-  PI ,)  <Xti  "W (156) 

tV-PO&3t+,W&«  (157) 


ajjxx  (v*  qw-  f?V+  QX,  +  PRi  a,,  +  t,V+ RV 

-PW-pi(+  Q.RA,')  a.,t  4-  [w+pv-  QU  (158) 
-lP,+  (Av)jl1]  a,. 
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(159) 


lu‘  +  QW-RV  +  <*-*-\+PR^)a.2l  +  (V4 -RU'-PW 

-PA,  +  QU,)aw  +  [w  +  pV-Q\T- 

0*^4,]  a„ 

QJ(Z*  ^GiW-RVt- GLl, +  PR.JL,)  ci3,  +(V+  RV  -PW 

-  Pi,  •I'QRL.V'iz  ■*■  [w+PV-QV-  (16°) 

(P**  <*')*'„]  «33 

The  numerical  solution  of  these  equations  then  represent 
the  desired  information  concerning  the  characteristics  of  the 
descent  trajectory  of  the  suspended  load  allowing  six  degrees 
of  freedom  of  the  parachute -load  system. 


2.  Summary  of  Assumptions  for  the  Free  Descent  Phase 

In  summary,  the  model  for  the  free  descent  phase  includes 
£he  following  assumptions.  (1)  The  parachute -load  system  can 
•be  described  as  a  rigid  body  whose  shape  does  not  change  as 
the  center  of  mass  location  changes  due  to  change  in  air 
density.  (2)  The  aerodynamics  of  the  parachute  can  be  fully 
described  by  a  tangent  force  and  by  normal  forces  and  aero¬ 
dynamic  moments  acting  at  the  moment  center.  (3)  The  aero¬ 
dynamics  of  the  load  are  described  sufficiently  by  the  drag 
force.  (4)  The  moments  about  the  center  of  mass  due  to  the 
suspension  lines,  risers,  extensions,  and  bridle  are  negligible, 
although  their  masses  are  included  in  calculation  of  the  center 
of  mass  position  and  the  moment  of  inertia  about  the  Y-axis. 

(5)  The  apparent  moment  of  inertia  can  be  sufficiently  approxi¬ 
mated  with  a  coefficient  measured  for  models  of  a  solid  flat 
circular  parachute,  and  the  apparent  mass  is  the  same  in  any 
direction. 


\ 
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IV.  NUMERICAL  TRAJECTORY  CALCULATIONS  OF 
SEVERAL  PARACHUTE-LOAD  SYSTEMS 


In  the  preceding  Sections  II  and  III,  various  mechanical 
arrangements  of  different  parachute  systems  and  their  function¬ 
ing  are  described  and  the  governing  equations  established. 

.In  view  of  the  overall  objective  of  this  study  a  continuous 
computer  program  incorporating  all  phases  of  the  functioning 
of  the  various  parachute-load  systems  was  developed.  This 
program  is  based  on  the  equations  shown  in  Section  III. 

In  order  to  check  the  functioning  of  this  program  and  to 
obtain  results  which  can  be  compared  to  a  certain  extent  with 
field  test  observations,  a  number  of  total  trajectory  simula¬ 
tions  were  carried  out.  The  results  of  these  calculations 
and  details  of  the  systems  for  which  the  simulations  were 
made  will  be  shown  later. 

At  this  time  it  may  be  stated  merely  that  the  program  . 
works  and  that  the  results  compare  favorably  with  field  test 
observations.  However,  before  reviewing  these  total  trajec¬ 
tory  simulations  it  appears  to  be  advisable  to  discuss  the 
more  general  inputs  to  the  governing  equations. 


A.  Steady  State  Canopy  Volume  and  Geometry 

Fcr  the  trajectory  simulations,  certain  assumptions  were 
made  concerning  the  shape  of  the  parachute  canopy  during  the 
free  descent  phase.  In  accordance  with  Ref  11,  the  volume  of 
a  fully  inflated  solid  flat  circular  parachute  is 

*7  s  (0.0GG7)  D03  (161) 

whereas  the  projected  diameter  of  a  fully  inflated  solid 
flat  circular  parachute  amounts  to 


°pwfcx  -  C°-6 °o 


(162) 


as  indicated  in  Ref  6.  Assuming  that  the  shape  of  the  para¬ 
chute  is  sami-ellipsoidal  and  the  volume  and  major  axis  are 
given  by  (161)  and  (162),  the  distance  to  the  parachute 
center  of  volume  from  the  canopy  skirt  amounts  to 

Sc.-(o.\Z9')  0o  (163) 

These  values  of  V,  and  sc  were  used  in  the  trajectory 

simulations  when  the  solid  flat  circular  parachutes  G-13, 

G-12D,  and  G-11A  were  used. 

For  the  extended  skirt  parachute,  type  T-10,  it  was 
assumed  that  (161)  is  valid,  whereas  measurements  of  photographs 
from  Ref  13,  combined  with  the  assumption  of  a  semi-ellipsoid/rL 
canopy  give  for  the  diameter  and  the  mass  center  distance 


(164) 

(165) 

The  geometry  of  the  parachute- load  suspension  system  is 
defined  by  the  value  of  the  lengths  of  the  suspension 

lines,  risers,  extensions  and  bridle,  and  the  load  dimensions. 
This  allows  the  determination  of  the  distances  s^  through  S5, 
and  together  with  sc  and  the  masses  of  the  components  of  each 
parachute- load  system  all  data  for  determination  of  the 
various  lengths  4  J[  2*  i  3*  an d  the  moment  of  inertia  lyy  are 
known,  provided  that  the  free  descent  phase  is  constrained  to 
three  degrees  of  freedom.  For  the  more  general  case  of  six 
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degrees  of  freedom,  the  only  additional  requirement  is 
sufficient  data  to  calculate  the  moment  of  inertia  about  the 
Z-axis,  l2z>  since  for  the  T-10  and  solid  flat  circular 
parachutes  the  moment  of  inertia  Ixx  *  lyy  311(4  Pr°duct 
of  inertia  Ixz  vanishes. 

Further  details  of  the  calculation  of  inputs  to  the 
trajectory  simulation  program  are  discussed  in  Volume  II  of 
this  report. 


B.  Aerodynamic  Coefficients 

The  remaining  information  necessary  for  the  trajectory 
simulation  are  the  aerodynamic  coefficients  for  the  force  and 
the  moment  for  the  T-10  and  solid  flat  circular  parachutes. 
Reference  14  presents  in  tables  and  graphs  the  results  of 
recently  conducted  measurements,  and  these  coefficients  were 
used  in  the  sample  calculations. 

The  tangent  force  coefficient  for  the  solid  flat  circular 
parachute,  taken  from  Ref  14,  can  be  approximated  by  the 
expression 

Ctq=  0.647 -(l.fcxlcT5’)  Kp|  +  pi 

-(.7.13  xio'OVXf!3  +(j-35xicft)l«f|1  O i  l°^p| *  Jo* 

’  '  1  (166) 

Ct,  =  O.CZO  ,  \ofy l  2  3o* 


Similarly,  the  tangent  force  coefficient  for  the  T-10  para¬ 
chute  is  approximately 

CTq=  0.570  ~(2.48xio'3)  Wp'  +(j-21')x|o"3)k(-|l 
-(7.4.87  x I o_r)  loCpl  5  +0.Z797  xio'^)  I«p|4  , 

OtUpl  C  (167^ 

Ct0  =  -.oo32(j<Xp|  -  30°K  .55-3  )  |c<pl230o 


The  measured  data  points  from  Ref  14,  and  these  approximate 
expressions  are  shown  in  Fig  15. 

The  aerodynamic  moment  coefficient  for  the  solid  flat 
circular  parachute  can  be  approximated  by 

Ca,q=  (4.<*44xio's)o(p  -(3.<H*IO-<1)o<p1L 

+  (,I.CH3  xle'Ootp  -d.34*lo‘‘,)o(p4  ,  OiotpOo'’ 

(168) 

CM.  =  -.0044-  Ctfp-'JO0  )-.03rt,  ot P  i  30° 


and  for  the  T-10  parachute  by 

*  (j.Wrf  x  \oz)  <Xp  -  xio'3>)  <*p 

-  (^8  xlo'5)  <*p3  -(.8.701) xio-7) o<(>4  (  Oioip^So° 

(169) 

cm.  =  “•oo6oUp-3e.*')-osfc,0<pi  30* 


Moment  coefficient  data  points,  Ref  14,  and  their  approxima¬ 
tions  are  shown  in  Fig  16. 

Finally,  the  normal  force  coefficients  for  the  solid 
flat  circular  parachute  can  be  ay proximated  by 

=  \0  °4p  1 1  >10 

-(i.s$xto'5)«P3  +(.;.9vio'7)o<p4  t  Oixfoo*(170) 

Ckj0  =  .ooSX  c  =<p-  3°*)  +.040,  <xp  >  3o* 
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SOLID  FLAT 
CIRCULAR 


Force  Coefficient  vs  Parachute 


(1.845  x  10-2)  _  (1.92  x  10-3)  a2 1  |  REF.  14  APPROX 


16  Moment  Coefficient  vs  Parachute  Angle  of  Attack, 
Solid  Flat  Circular  and  T-10  Parachutes 


and  for  the  T-10  parachute  by 

Cw&=  -(2.5'o8  xlO~2)*p  + O-^'Cf^Vp 

-(^>•02^K\0"S■,)  «*  (jo.8  2*7  X  \67)  <tf?4  ,  o  <  0(p<  3o° 

(171) 

CWft*  .0072.  Co(p-  3o*)  +  .(xA;<*p  £  30  * 


Again,  data  points  from  Ref  14  and  expressions  (170)  and  (171) 
are  shown  in  Fig  17,  Because  of  the  parachute  symmetry,  the 
signs  of  the  coefficients  of  C^0  and  Cj^  are  reversed  for 
negative  angles  of  attack. 

It  should  be  noted  that  the  signs  for  the  normal  forces, 
t  Fig  12,  and  coefficients,  Fig  16,  are  reversed  from  those 

,  shown  in  Refs  14  and  15  in  order  to  arrange  them  in  a  standard 

right-hand  coordinate  system.  This  was  done  because  it  was 
|  necessary  that  in  this  trajectory  simulation,  all  forces, 

moments,  and  angles  follow  a  sign  convention  consistent  with 
a  right-hand  coordinate  system. 


r 


C.  Examples  of  Total  Trajectory  Simulations 

Once  the  governing  equations  had  been  formulated  and  the 
computer  programming  completed,  several  sets  of  input  data 
were  selected  to  provide  a  basis  for  debugging  the  program, 
demonstrating  that  all  of  the  calculation  methods  function  as 
intended,  and,  of  course,  to  see  if  the  results  agreed  reason¬ 
ably  well  with  field  test  experience.  Five  parachute -load 
systems  were  used  with  various  release  velocities,  release 
altitudes,  suspended  loads,  and  main  parachute  reefings  to 
give  inputs  for  a  data  production  run  of  21  total  trajectory 
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ig  17  Norma!  Force  Coefficient  vs  Parachute  Angie 
Attack; Solid  Flat  Circular  and  T-'IO  Parachutes 


simulations  specified  by  the  procuring  agency.  The  input 
data  for  these  21  cases  is  presented  in  detail  in  the 
Appendix.  The  data  was  not  intended  to  provide  the  extensive 
parametric  variations  which  would  be  needed  to  judge  the 
accuracy  of  the  calculations. 

All  of  the  sample  calculations  presented  were  constrained 
to  three  degrees  of  freedom,  since  no  measurements  of  aero¬ 
dynamic  coefficients  for  a  system  allowing  six  degrees  of 
freedom  exist.  However,  one  sample  calculation  was  accom¬ 
plished  with  six  degrees  of  freedom  using  the  same  aero¬ 
dynamic  coefficients  in  the  XZ-  and  YZ-planes.  These  functions 
were  the  same  as  those  used  with  three  degrees  of  freedom. 

Since  no  disturbance  in  the  YZ-plane  was  introduced,  the 
calculation  gave  results  which  remained  two-dimensional  and 
were  exactly  the  same  as  for  the  calculation  allowing  only 
three  degrees  of  freedom.  This  indicates  that  the  total 
trajectory  simulation  with  six  degrees  of  freedom  functions 
properly. 

In  the  following,  graphical  and  numerical  results  are 
presented  for  calculations  for  four  separation -deployment 
methods,  including  ore  case  with  disreefing  of  the  main 
parachute.  Also,  some  results  for  the  remainder  of  the  21 
sample  calculations  are  tabulated. 


1.  T-10  Parachute  with  Static  Line  System  • 

The  physical  data  for  this  calculation  were  from  a  T-10 
personnel  parachute,  D0  =  35  ft,  with  a  suspended  load  of 
250  lb.  In  this  case  it  was  released  from  an  altitude  of 
6,000  ft  at  a  speed  of  220  ft/ sec.  The  physical  processes 
which  occur  throughout  the  airdrop  for  the  static  line  system, 
followed  by  the  particular  Subsection  or  Part  of  this 
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report  in  which  the  modeling  equations  were  developed,  are: 
load  falling  away  from  aircraft,  IIIC2,  deployment  of  main 
parachute,  IIID1,  main  parachute  inflation,  I  HE,  and  free 
descent,  IIIF. 

The  trajectory  variables  of  altitude  loss,  z,  system 
angle,  0,  horizontal  velocity,  and  vertical  velocity  are 
plotted  versus  time  in  Figs  18  and  19.  These  show  a  system 
oscillation  with  a  period  of  about  6  sec,  which  has  to  be 
expected  since  the  T-10  is  unstable  at  =  0°.  The  small 
discontinuity  in  the  curves  at  3.2  sec  is  the  result  of 
changing  from  data  at  the  mass  center  to  data  at  the  load 
center  after  full  inflation.  This  discontinuity  occurs  in  all 
cases,  and  becomes  larger  as  the  parachute-load  system  lengths 
increase.  For  this  T-10  static  line  system  there  is  no  snatch 
force  calculation,  the  unreefed  filling  time  was  1.65  sec, 
and  the  maximum  opening  force  was  1868  lb. 


2.  G-12D  Cargo  Parachute  with  Static  Line  Deployed 

Pilot  Chute  System. 

The  parachute- load  system  that  provided  the  physical 
data  for  this  calculation  was  a  G-12D  parachute  with  a  sus¬ 
pended  load  of  2200  lb.  The  G-12D  is  a  solid  flat  circular 
parachute,  DQ  ~  64  ft,  also  unstable  at  Op  -  0°.  The  release 
conditions  were  again  hQ  13  6,000  ft,  vq  *  220  ft/sec.  The 
physical  processes  and  modeling  equations  were  :  load  falling 
away  from  aircraft,  IIIC2,  suspension  system  deployment, 
IIID2a,  snatch  force,  IIID2b,  main  parachute  unfolding, 

IIID2c,  main  parachute  inflation,  IIIE,  free  descent,  IIIF. 

The  trajectory  variables  calculated  for  this  case  are 
shown  in  Figs  20  and  21.  Again  a  damped  oscillation  is  cal¬ 
culated,  but  with  a  system  larger  than  the  T-10  system,  the 
period  is  about  8  sec.  The  figures  show  that  this  system  is 
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Altitude  Loss  and  System  Angle  for  the  T-10 
Parachute  with  Static  Line  System 


ho  =  6000  ft 


Fig  20  Altitude  Loss  and  System  Angle  for  the  G-12D 
Cargo  Parachute  with  Static  Line  Deployed 

Pilot  Chute  System 


=  6000  ft 
=  220  (ft/sec) 

=2200  lb 


Deployed  Pilot  Chute  System 


approaching  a  steady  glide  related  to  the  parachute's  stable 
angle  of  attack,  as  all  systems  did  when  given  enough  alti¬ 
tude.  The  snatch  force  was  8293  lb  at  a  velocity  of  197  ft/ 
sec,  the  unreefed  filling  time  was  1.65  secs,  with  a  maximum 
opening  force  of  11,847  lb. 


3.  G-11A  Cargo  Parachute  with  Extraction  Parachute  System. 

a.  Unreefed  Main  Parachute  Inflation 

Data  for  this  case  was  obtained  from  a  G-11A  parachute, 
Dq  =  100  ft  with  a  3500  lb  suspended  load.  The  extraction 
parachute  used  was  a  15  ft  reefed  ringslot  parachute.  The 
release  conditions  were  hQ  =  2,000  ft,  vQ  *  220  ft/ sec.  The 
physical  processes  and  modeling  equations  were:  extraction 
of  load  from  aircraft,  IIIC1,  suspension  system  deployment, 
IIID2a,  snatch  force,  IIID2b,  main  parachute  unfolding, 
IIID2c,  main  parachute  inflation,  HIE,  and  free  descent, 
IIIF. 

The  trajectory  variables.  Figs  22  and  23,  again  show  a 
damped  oscillation  with  a  period  of  11  secs.  The  calculated 
snatch  force  was  2,357  lb  at  a  velocity  of  160  ft/ sec,  and 
the  maximum  opening  force  was  9,331  lb  for  an  unreefed  fill¬ 
ing  time  of  5.43,  sec. 

b.  Reefed  Main  Parachute  Inflation 

This  calculation  was  made  to  demonstrate  that  reefed 
inflations  could  be  modeled.  The  case  presented  here  uses 
the  same  separation- deployment  system  and  parachute- load 
system  as  discussed  in  the  previous  item.  The  release 
conditions  were  hQ  =  2,000  ft,  vQ  =  169  ft/ sec.  This  sample 
case  represents  a  three  stage  reefed  inflation  with  reefing 
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G-11A  Cargo  Parachute  with  Extraction  Parachute 
System 
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Fig  23  Horizontal  and  Vertical  Velocities  for  the 
Unreefed  G-11A  Cargo  Parachute  with 
Extraction  Parachute  System 


cutter  delay  times  of  2  sec,  4  sec,  and  6  sec  after  snatch 
force.  The  equations  of  subsection  IIIB  were  used  for  the 
coasting  periods  of  constant  drag  area  between  inflation 
stages. 

Figures  24  and  25  present  the  trajectory  variables,  and 
again  show  a  damped  oscillation  with  an  11  sec  period.  It 
is  interesting  that  the  amplitude  of  the  oscillations  was 
considerably  smaller  than  the  unreefed  G-11A.  The  snatch 
force  was  2078  lb  at  a  velocity  of  130  ft/sec;  the  reefed 
inflation  times,  excluding  cutter  delays,  were  1.68,  0.96,  0.72, 
and  4.01  sec  with  corresponding  maximum  forces  of  3,062  lb, 

5,667  lb,  6,032  lb,  and  5,828  lb. 

4.  G-11A  Cargo  Parachute  with  Reefed  Main  Parachute 

Extraction  System. 

The  physical  data  for  this  system  came  from  a  G-11A 
parachute  with  a  20- ft  reefing  line  and  a  3500  lb  extracted 
and  suspended  load.  The  reefing  line  was  cut  4  sec  after  the 
load  was  released  in  the  aircraft.  The  release  conditions 
were  hQ  ■  2,000  ft,  vQ  *  220  ft/ sec.  This  sample  case  has 
the  following  physical  processes  and  modeling  equations: 
extraction  from  aircraft,  IIIC1,  main  parachute  inflation, 

HIE,  and  free  descent  IIIF. 

The  trajectory  variables  are  presented  in  Figs  26  and  27 
and  were  very  similar  to  the  other  G-11A  systems.  The  maxi¬ 
mum  opening  force  was  8,348  lb  for  a  reefed  inflation  time 
of  5.12  sec. 

5.  Additional  Calculations. 

Tables  I  and  II  show  numerical  results  for  all  of  the 
21  total  trajectory  calculations  made,  including  those  just 
discussed.  Graphs  like  those  presented  previously  were  made 
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Altitude  Loss  and  System  Angle  for  the  Reefed 
G-11  A  Cargo  Parachute  with  Extraction  Parachute 
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Fig  25  Horizontal  and  Vertical  Velocities  for  the 
Reefed  G-11A  Cargo  Parachute  with 
Extraction  Parachute  System 


2000  ft 


Fig  26  Altitude  Loss  and  System  Angle  for  the  G-11A 
Cargo  Parachute  with  Reefed  Main  Parachute 
Extraction  System 


Parachute  Extract  ion  System 


H 

lli 

s 

2 


PARACHUTE -LOAD  SYSTEMS  AND  INITIAL  CONDITIONS 

FOR  THE  TOTAL  TRAJECTORY  SIMULATION  CALCULATIONS 

V 

>«* 

b 

W 

o 

N 

N 

o 

in 

n 

■* 

*4 

o 

N 

fN 

V7 

m 

+4 

8  8 

N  «* 

o  a 

N  A 
N  « 

o 

N 

IN 

o 

IN 

fN 

r— 

o 

IN 

IN 

O 

fN 

IN 

o 

IN 

IN 

o 

IN 

IN 

o 

IN 

IN 

o 

'N 

fN 

O 

fN 

IN 

O' 

40 

>4 

o 

IN 

IN 

o 

IN 

IN 

O  *? 
£  <M 

w 

s 

V 

o 

o 

ia 

S 

O 

8 

i  i 
« 

o  a 

8  S 

IN 

8 

o 

N 

s 

Q 

IN 

8 

c 

o 

r4 

O 

O 

O 

IN 

O 

o 

o 

o* 

r4 

o 

8# 

IN 

o 

o 

o 

4 

o 

-4 

8 

c 

IN* 

o 

r4 

8 

o 

IN 

IN 

O 

o 

o 

IN 

0  « 

*  : 

1 

i 

1 

1 

• 

• 

i 

< 

IN 

IN 

NT 

IN 

IN 

* 

« 

NT 

-■> 

IN 

4f 

sp 

m  tf 
**c  a, 

« 

i 

« 

1 

< 

i 

4 

< 

8 

8 

o 

IN 

8’ 

8 

O 

IN 

O 

fN 

8 

8 

"s 

S 

O 

IN 

O 

O 

IN 

** 

.-a 

8 

N 

8 

N 

O 

8 

« 

S 

lA 

^  f 

N  f> 

8  S 

IN  W 
*  IN*  P 

’  A 

8 

lA 

n" 

8 

1" 

O 

O 

tA 

n 

<a 

8 

o 

lA 

O 

8 

WN 

o 

8 

lA 

•A 

8 

<A 

fl* 

IA 

i 

n* 

l 

85 

I 

s 

•4 

V 

ik 

3 

M 

! 

a 

V 
*4 
u 
• 
u 

VI 

S 

*4 

~| 

a 

u 

3 

w 

! 

-4 

>4 

V 

■4 

b 

• 

b 

M 

s 

•J 

V 

b 

W 

• 

u 

M 

• 

5  f 

b 

3  . 

£  1 

1 

H  0 

a 

2 

2 

-j 

V 
•4 
b 
• 

W 

VI 

1 

4 

4 

4 

• 

44 

2  ■* 

i  • 

z. 

s 

u 

Si  i 

4  b 
b 

X 

*1  N 

5  3 

1  1 

£ 

£ 

5  8 

2  t 

i  2 

k  14 

-»  w 

3  £ 

• 

44 

a 

2 

b 

£ 

8 

vH 

u 

2 

b 

w 

M 

M 

■ 

b 

£ 

I 

8 

b 

b 

U 

• 

b 

44 

X 

M 

Cl 

4k 

£ 

2 

b 

£ 

i 

k 

4k 

U 

c 

V 

4k 

K 

w 

2 

1 

s 

£ 

8 

k 

2 
b 
b 

X 

M 

O 

b 

1 

2 

■ 

lb 

s 

—4 

44 

U 

■ 

b 

b 

£ 

2 

a 

2 

•  b 

£ 

8 

k 

4k 

u 

* 

b 

u 

£ 

2 

£ 

u 

b 

£ 

§ 

—i 

b 

u 

• 

b 

44 

X 

M 

2 

a 

i* 

2 

§ 

-4 

b 

u 

■ 

b 

b 

£ 

• 

b 

2 

1 

£ 

8 

-4 

b 

a 

u 

b 

b 

X 

M 

8 

4k 

u 

« 

k 

4k 

£ 

2 

£ 

2 

b 

£ 

c 

•4 

£ 

•  • 

M 

3 

£ 

2 

b 

y 

« 

M 

b 

K 

U 

2 

£ 

1 

c 

•vk 

2 

• 

N4 

3 

! 

o 

r* 

|4 

s 

0 

lA 

r4 

4 

lA 

r4 

A 

« 

0k 

A 

a 

IN 

44 

6 

a  e 

4  IN 

■4  ~i 

i  o 

<  < 
■4  *4 

4  -4 

5  o 

«C 

0k 

0k 

t!i 

< 

r4 

*4 

A 

< 

pk 

+4 

A 

< 

0k 

4 

< 

0k 

» 

o 

< 

0k 

» 

o 

< 

r4 

-4 

• 

u 

< 

0k 

0k 

4 

< 

0k 

0k 

u 

a 

0k 

A 

a 

r4 

6 

i 

A 

o 

«* 

>A 

o  ! 

t  - 

!»  o 

*4 

|4 

r4 

IN 

r* 

ia 

44 

if 

1-4 

lA 

r4 

N» 

*4 

A. 

44 

CO 

0k 

S' 

0k 

& 

fN 

r4 

IN 

83 


*Graphical  Results  Presented 


for  all  cases,  but  they  have  not  been  reproduced  because 
they  were  essentially  the  same  as  those  presented.  They 
showed  differences  in  oscillation  period  and  terminal 
velocity,  but  this  would  be  expected  from  systems  with  various 
sizes  and  loadings.  The  accuracy  of  the  results  cannot  be 
judged  at  this  time,  but  obtaining  all  of  these  results  from 
a  single  submission  of  the  computer  program  certainly  shows 
that  the  total  trajectory  simulation  functions  satisfactorily 
with  any  pertinent  inputs. 


V.  SUMMARY 


The  functioning  of  an  airdrop  or  recovery  system  by 
me^os  of  parachutes  was  organized  in  phases  covering  the 
main  events.  These  were  the  separation  of  the  load  from  an 
aircraft,  the  deployment  of  the  parachute  system,  the  dynamic 
and  aerodynamic  forces  during  the  parachute  snatch  and  para¬ 
chute  inflation  processes,  and  the  movement  of  the  parachute - 
load  system  after  the  inflation  of  the  main  parachute.  The 
phases  indicated  above  included  considerations  related  to  the 
various  methods  required  by  the  differing  arrangements  of  the 
parachute  systems.  For  example,  consideration  was  given  to 
the  inflation  of  the  main  parachute  in  an  unrestrained  single 
process  as  well  as  in  steps  accomplished  by  different  reefing 
line  arrangement?. 

For  all  phases,  the  basic  equations  were  taken  from 
existing  literature,  when  available,  or  they  were  established, 
as,  for  example,  the  equations  of  motion  of  the  fully  in¬ 
flated  and  descending  parachute  allowing  3ix  degrees  of 
freedom. 

A  continuous  computer  program  was  developed  which  allowed 
the  total  trajectory  simulation  of  any  parachute -load  system 
whose  arrangement  fits  any  one  of  the  cases  discussed  in  this 
study.  The  program’,  was  designed  in  such  a  way  that  different 
subroutines  relating  to  mechanical,  aerodynamic,  and  dynamic 
methods  can  be  introduced,  removed,  or  exchanged. 

After  completion  of  the  computer  program  a  number  of 
numerical  total  trajectory  calculations  was  carried  out  and 
the  results  were  indicated  in  tables  and  graphical  presenta¬ 
tions.  A  study  of  the  calculated  results  showed  that  they 
fell  well  within  the  broad  ranges  one  would  expect  based  upon 
a  familiarity  with  existing  field  test  results.  All 
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calculated  results  follow  the  known  basic  parachute  performance 
characteristics  very  well.  Statements  about  the  accuracy  of 
the  calculated  results  cannot  be  made  without  having  more 
complete  field  test  data.  Also,  the  calculation  of  the  move¬ 
ment  of  the  parachute- load  system  after  full  inflation  of  the 
main  parachute  had  to  be  restricted  to  movement  in  one  plane 
because  aerodynamic  coefficients  of  parachutes  at  composite 
angles  of  attack  were  not  available. 

The  computer  program  has  been  laid  out  and  calculations 
have  been  made  which  consider  the  known  non-linearities  of 
the  aerodynamic  coefficients  versus  their  angle  of  attack. 

The  calculations  do  prove  that  the  complete  single  computer 
program  is  capable  of  predicting  the  dynamics  of  the  entire 
aerial  delivery  process.  In  particular  it  has  been  shown 
that  the  movements  of  the  parachute- load  system  in  two 
orthogonal  planes  can  be  predicted  provided  that  the  aero¬ 
dynamic  inputs  are  available. 

It  appears  that  the  program  is  ready  to  be  used  for 
predicting  parachute  performance  and  that  intensive  compari¬ 
sons  between  calculated  and  observed  performance  character¬ 
istics  are  highly  desirable. 
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APPENDIX 


The  numerical  data  for  the  21  sample  trajectory  cal¬ 
culations  discussed  in  Section  IV  were  based  on  five  parachute¬ 
load  systems  frequently  used  by  the  United  States  Army.  The 
parachutes  used  in  these  five  systems  are  the  T-10  personnel 
parachute  and  the  G-13,  G-12D,  and  d-llA  cargo  parachutes. 

The  four  separation -deployment  systems  discussed  in  Section 
II  correspond  to  these  systems  as  follows:  the  static  line 
system  (System  1)  is  used  for  the  T“10  and  G-13  parachutes, 
the  static  line  deployed  pilot  chute  system  (System  2)  is 
used  for  the  G-12D  parachute,  the  extraction  parachute  system 
(System  3)  is  used  with  the  G-11A  and  G-12D  parachutes,  and 
the  reefed  main  parachute  extraction  system  (System  4)  is 
used  with  the  G-11A  parachute.  All  physical  details  of  these 
five  parachute- load  systems  as  used  with  the  four  separation- 
deployment  methods  are  indicated  schematically  and  described 
in  Figs  28  through  32.  The  various  masses,  lengths,  volumes, 
etc.  which  were  used  as  inputs  to  the  total  trajectory 
simulation  were  determined  from  Figs  28  through  32  and  are 
listed  in  Tables  Kl  through  VII.  The  drag  coefficient  Cj)p  is 
given  by 

v/  •' '  • 

f-'  ^  Oo 

i-Dp= - 

(DPwy  /Oo)Z 
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1)  lode  Um,  iwtU  •  13  ft,  MUriil  •  XXLHMOM,  Iff  VtU 
•ylaa. 

Z )  BipUjMU  k|,  laagefc  •  22  la,  vlitk  •  12  la,  balfht  • 

2  la,  (Cgf)  «*  0.33  ft1,  «tl«hc  •  3.3  Ik. 

3)  Braak  aanl,  aatarlal  -  Typa  1,  k  wMlit,  cattaa,  finish^  ’aagtfc 
•  3  la,  vvtffee  •  0.11  m. 

4)  *#  •  S3  ft,  10X  flat  uttoM  aklrt  «aww.  Caaapy  aatarlal, 
aylau  alath  KIL-C-7020,  Tyya  1.  1.1  aa/jNT,  aaalaal  poroalty  -  72  - 

132  ft’/f^-aia  (um  100  ft*/f  t2-«U) .  W.  -  11.7  lk,  I  - 
14.13  lk,  C-  -  0.70,  *  * 

% 

3)  30  saaaaaflaa  llaaa,  KIL-C-5040,  Tjrpa  IX,  373  lk,  •  32X, 

1(  -  23^3  ft,  wl|bt  •  2.42  lk.  “  . 

4)  4  rtaara.  aatarUl  -  XIL-*-«*4,  Tjrpa  ZUI  aylaa,  -  301  (aat), 

•  30  la,  walfht  •  1.0  lk, 

7)  toylaaf  V,  •  230  lk,  (Cgl),  •  4.0  ft1. 


Fig  28,  T-10  Parachute  with  Static  Line  System 


TABLE  III 


TOTAL  TRAJECTORY  SIMULATION  INPUTS 
FOR  T-10  PARACHUTE  WITH  STATIC  LINE  SYSTEM 

Dimensions  and  Lengths 

D  “  35  ft  D  -  (0.6S6)D  sn  -  11.52  ft 

o  n  o  i 

nnax 

Ls  -  25  ft  sc  -  (0.133) Dq  s2  -  24.17  ft 

^R  *  2-5  ft  ^static  "  15  ft  8 5  "  27‘8  ft 
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1)  Stacie  Una,  Static  "  **  *c  lo»  “  MIL-W-5665,  “type 

VIII  cotton  web,  2900  lb  or.  ftr.  sin. 

2)  Braak  cord,  caterial  “  Type  I  k  Inch  cottcn  vebblog,  80  lb, finished 
length  “  6  In,  weight  "  0.17  ox. 

3)  DQ  •  32  ft  heu l spherical  canopy,  canopy  Material,  cotton  kusIIo, 

4.25  o z/yd2  KIL-C-4279,  Type  il,  ooolnal  poroalty  •  153  -  253  ft*/ 

ft2-tdo,  V  -  23. £9.  V  •*  33.69,  C„  -  0.75. 

p  o 

4)  20  suspension  lines,  MIL- C* 4 2 32  Type  I,  40C  lb,  ■  12%,  * 

30  It,  weight  -  4.8  lb. 

5)  Risers,  tutorial  •  KIL-V-S665.  Tyo*  VIU  cotton  vsfcblng  2900  lb  (aio), 

€  mAX  ■  8  -  ICZ  (ast),  Lg  -  30  lo,  weight#*  1.2  lb. 

6)  Cargo  container,  length  -  30  In,  width  •  X  in,  height  -  20  In, 
type  A-21,  (Cp*),  -  7.125  ft2,  W4  -  500  lb. 


Fig  29.  G-13  Cargo  Parachute  with  Static  Line  System 
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TABLE  IV 


TOTAL  TRAJECTORY  SIMULATION  INPUTS  FOR  G-13 
CARGO  PARACHUTE  WITH  STATIC  LINE  SYSTEM 


Dimensions  and  Lengths 


Dq  =  32  ft 

D  -  (0 . 648)  D 

‘max 

sx  -  14.215  ft 

Ls  =  30  ft 

sc  “  (0*129) Dq 

s2  -  29.615  ft 

1^  -  2.5  ft 

^static  =  15  ft 

s5  «  31.635  ft 

Masses 

nip  ■  0.743  slug 

*  0.037  slug 

m  -  16.471  slug 

i  s 

nw  *  0.149  slug 

Ls 

mj  -  15.542  slug 

Drag  Area  and  Drag  Coefficient 

CnS.  -  7.125  ft2  CD  -  1.786 

%  P 


Main  Parachute  Canopy  Volume 
V  -  2185.62  ft3 
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1) 

u 

3) 

4) 

3) 

4) 

7) 

V. 

9) 

M> 

ID 

It) 


C-12D  Cargo  Parse  hue  a 


Stacie  11m,  Lit<tic  “  15  ft,  aatarlal  “  Type  VIZZ  cotton  webbing, 
rolled  KXL-W-5665. 

Pilot  PnacHuta  Pack,  might  •  0.5  lb,  (CgfJwO.S  ft2. 

Pilot  Paroehuta.  D„  "  5.66  ft,  octagonal  canopy.  3  auaoaroalon 
llBM.  L  -  5. 5  ft?  KZL-C-7S1J,  Typ*  ZZ,  C~  •  0.75,  valght  - 
1,75  lb,  o 


Deeloynant  lloa  and  brldla,  notarial  “  Typo  VIII  nylon  webbing 
^■ax  *  (••*)»  ltngch  -  14  ft,  walght  **  2.0  lb  w/clevla. 


Deplcyvart  bag,  dinenatooo,  laogth  •  36  lo,  width  “  24.5  lo,  height  “ 
10  lo,  (CpS)  -  1.78  ft2,  might  -  6  lb. 

Braak  cord,  aatarlal  •  Type  IZZ  550  lb  nylon  cord  MIL-C-5040, 
flotahed  laogth  -  6  In,  weight  “  0.005  lb. 


•  64  ft  .c-rcular  flat  canopy,  canopy  material,  nyloo,  2.25  ns/ 

KIL-C-7350,  type  I.  float nil  pcroalty  -  50  -  165  ft3/ft2-«ln, 

•  74.79  lb,  -  1C§.5  lb. 


64  auaoenaloo  llnea,  KH.-C-7315,  Type  IV,  1,000  lb.  -  20X, 
l,  -  51  ft,  weight  -  27.9  lb.  *■* 

f  r  leer  a,  KIL-W-27265,  Type  X,  Claaa  t,  -  25X,  8700  lb, 

4  -  5  ft  5  in,  weight  -  4. 7? lb.  “* 


4  Sutpetialoo  SUoga,  KXL-V-27265  Type  VZZ  oyloo  web,  5500  lb, 
•  »X,  length  -  26  In,  weight  «2  lb. 


16  A-22  Suapenaton  Sllnca,  KZL-W-5665  Tyne  VZlt  cotton  veb,  2900  lb, 
*ivax  *  *****  lfn8ch  *  5  ft,  weight  7  lb  (ahown  on  front  aide  of 
container;  typical  of  all  four  aides) • 


A  22  Cargo  coataloar,  length  *  52  Ln,  width  -  43  to,  height  "  50  In, 
(CB*)J  -  18.9  ft2,  Vj  -  2,200  lb  (with  cargo) 


Fig  30.  G-12D  Cargo  Parachute  with  Static  Line 

Deployed  Pilot  Chute  System 
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TABLE  V 

TOTAL  TRAJECTORY  SIMULATION  INPUTS  FOR  G-12D 
CARGO  PARACHUTE  WITH  STATIC  LINE  DEPLOYED 
PILOT  CHUTE  SYSTEM 


D  =  54  ft 
o 

Lg  =  51  ft 
1^  =  5.42  ft 

Ljir  “  8  ft 


Dimensions  and  Lengths 


Dd  *  (0.648)D 
Fmax 

•c  “  <°-129>Do 


D 


'static 
o) pilot 


<Ls  +  pilot 


15  ft 
5.66  ft 
14.75  ft 


a. 


23.75  ft: 
50.00  ft 
56.23  ft 
57.21  ft 


m^  -  2.383  slug 
°  0.867  slug 
mR  -  0.147  slug 


Masses 

mBr  -  0.280  slug 
m ^  *  68.386  slug 


pb 


Drag  Areas  and  Drag  Coefficient 


CDS.  -  18.9  ft2  CDSn  -  26.016  ft2  CD 

*  P 


0.303  slug 
72.870  slug 


1.786 


Main  Parachute  Canopy  Volume 
V  -  17,301.5  ft3 

Suspension  System  Spring  Constant 
1/A  -  4777  lb/ft 
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C-120  Cargo  Parachute 


l)  15  ft  ringiloe  parachuc*  -  reefed  Co  D 


2) 


eff 


-  12  fc,  caoopy 

,2 


material,  nylon  cloth  MIL- C- 7 350 ,  Type  I,  2.25  oz/yd  ,  noaical 
poroaity  -  90  -  165  ft3/ft2-*in,  Cj,  •  0.55,  Wc  -  6.06  lb,  Wf 
M6  lb,  <V>r,.,td  -  «1.1  *tJ.  0 
Reefing  lint,  lcngch  ■  260  in. 


3)  Extinction  parachuc*  brldla;  materiel  -  WL-W-27265.  nylon,  Typt  X, 
1700  lb  br.  at.  (min),  C  ■  301  (ate),  length  -  60  fc,l  loop, 
Wliht  »  10  lb.  *“ 


4) 

5) 

6) 

7) 

*) 


Daploycanc  bag,  dimension*,  length  “  36  In,  vidcb  •  24.5,  halghc  " 
10  lo,  (CpS)  -  1.73  ft2,  weight  -  6  lb. 

IreaV  cord,  material  -  Tyne  III  550  lb  nylon  cord  KH.-C-5040, 
fioi:h*d  langch  •  6  in,  weight  •  0.005  lb. 


0O  -  64  f t  circular  flat  canopy,  canopy  material,  nylon,  2.25  ox/ 

Jd2  KlL-C-7350,  Type  I,  nemtoal  poroiity  -  90  -  165  ft3/fe2-mio, 

#  -  76.79  lb,  Wp  -  10$. 5  lb. 

64  auapeniioo  linn,  KtL>C>75Lr,,  Type  IV,  1,000  lb,  (  -  20X, 

l#  -  51  fc,  weight  -  27.9  lb. 

8  riaere,  KIL-V-27265,  Tyne  X,  Clai*  t,  -  28X,  8700  lb,  L,  • 
5  ft  5  in,  weight  -  4.75  lb.  * 


9)  fUtfora  bridle,  material  -  MIL-W-27265,  nylon  Tabbing,  Type  X, 
8700  lb  br.  etr.  (tin),  CMX  *  30X  (eat),  length  -  9  fc  (4  each 

2  loop)  weight  -  11,0  lba. 

10)  Cargo  platform,  length  -  8  ft,  width  ■  8  ft,  (CqS)*  •  76.3  fc2, 

W,  *  2,200  lb  (with  cargo) •  9 


Fig  31.  G-12D  Cargo  Parachute  with  Extraction 

Parachute  System 

% 
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TABLE  VI 

TOTAL  TRAJECTORY  SIMULATION  INPUTS  FOR  G-12D  CARGO 
PARACHUTE  WITH  EXTRACTION  PARACHUTE  SYSTEM 

Dimensions  and  Lengths 

Lfir  -  8  ft  s1  -  23.75  ft 

D  =  (0.648)D  s2  =  50.00  ft 

^max  0  z 

s  -  (0.129)D  s,  =  56.50  ft 

C  O  4 

s5  *  62.50  ft 

Masses 

mp  =  2.383  slug  m^  -  0,341  slug  mpb  *  0.746  slug 

m^  *  0.867  slug  m^  *  68.386  slug  mrg  =  72. 

s 

mR  -  0.147  slug 

Drag  Areas  and  Drag  Coefficient 
CDSjt=  76,8  ffc2  CDSII  “  63,98  ft2  CD  “  1'786 

Main  Parachute  Canopy  Volume 
V  -  17,301.5  ft3 


D  =  64  ft 
o 

Ls  =  51  ft 
L^  ■  5.42  ft 


Suspension  System  Spring  Constant 
1/A  -  4777  lb/ft 


1) 


2) 

*) 

*) 

3) 

«) 


*) 


IS  ft  rirgilot  parachute  -  raafad  to  D  -  12.0  ft,  canopy 

°eff 


material,  nylon  cloth  HI1-C-7350,  Typa  I,  2.2S  o*/yd  ,  non  Inal 
poroalty  -  90  -  165  ft3/ft2-mln,  CD  -  0.55,  Vc  -  6.06  lb,  II 
*.06  lo.  <C0Sl  ,  .  -  62.2  ft2.  0 


'O'*' reefed 
leafing  J.loa,  length  -  260  la; 


Extraction  parachute  brldl*;  material  -  KIL-W-27265  nylon,  Typa  X, 
*700  lb  br.  att.  (dn),  €  -  301  (aat),  langth  -  60  ft  (one 

loop),  weight  •  10  Iba. 

iVployaan t  bag,  dioenaions,  langth  *  48  In,  width  -  35  In,  height  • 
12  In,  (CDS)  -  2.33  ft2,  weight  •  28  lb.. 

Break  cord,  oatar ial  -  Typa  IH  550  lb  nylon  cord  KT1.-C-5040, 
finished  langth  «  6  in,  weight  •  0.005  lb. 


0o  -  100  ft  circular  flat  canopy,  canopy  notarial,  nylon,  1.6  ox/ 

yd2  -MIL- C- 7020.  Type  II,  nominal  poroalty  *•  70  -  176  ftVft2-*i«, 
Wf  -  133.7  lb,  Wp  -  193  Lb. 


120  suspension  lines,  KIL-C-5040,  type  III,  550  lb,  r  351, 
L,  -  35  ft,  weight  -  20.7  lb.  *  “•* 


12  risers,  MIL-W-27265,  Type  XVIII,  6,000  lb,  -  211,  L»  - 
60  It,  weight  -  33.6  lb.  ®*x  ^ 


9)  Load  attachment  riae.,  material  -  MlL-W-27265  nylon,  Typa  X, 

*700  lu  br.  str.  (oinj,  (  -  301  (*«t),  L.  “  20  It  (one  loop), 

weight  -  4  lb.  ***  6 


10) 

ID 


Platform  bridle,  rcaterial  •  KIL'W-27265  nylon.  Type  X,  8700  lb 
br.  atr.  (min),  ■  XI  (ast),  length  -  9  ft’,  weight-  11. 0 

Cariio  plat  fora,  length  -  5  ft,  width  -  6  ft,  (CnS).  -  76.8  ft2, 
V.  -  3500  lb  or  500u  lb  (with  cargo) <  0 


lb. 


Fig  32.  G-11A  Cargo  Parachute  with  Extraction 

Parachute  System 
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TABLE  VII 


TOTAL  TRAJECTORY  SIMULATION  INPUTS  FOR  G-11A  CARGO 
PARACHUTE  WITH  EXTRACTION  PARACHUTE  SYSTEM  OR 
REEFED  MAIN  PARACHUTE  EXTRACTION  SYSTEM 

Dimensions  and  Lengths 

DQ  =  100  ft  Lfir  -  8, ft  s±  «  16.‘45  ft 

L  -  35  ft  D  -  (0.648)D  -  61.10  ft 

s  pmax  °  4 

1^-  60  ft  sc  -  (0.129)Dq  s 3  -  99.30  ft 

s4  -  113.30  ft 

s5  -  119.30  ft 

Masses 

mp  -  4.156  slug  Qgr  -  0.342  3lug  mpb  *  1.432  slug 

nw  ■  0.643  slug  m§  ■  108.8  slug  m  *  116.697  slug 

s  *  or  rs  or 

155.4  slug  163.324  slug 

mE  -  0.124  slug 

Drag  Areas  and  Drag  Coefficient 
CDSj  -  76.8  ft2  CDSn  -  64.53  ft2  Cj,  -  1.786 

Main  Parachute  Canopy  Volume 
V  -  66,700  ft3 

Suspension  System  Spring  Constant 
1/A  -  910.7  lb/ ft 
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